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The electric field in mm-sized clouds of cold Rb + -ions is 
measured using the Stark effect of Rydberg atoms embed- 
ded in the ion clouds. The ion clouds are produced by UV 
photoionization of laser-cooled rubidium atoms in a magneto- 
optic trap. In very dense ion clouds, the Rydberg atom ex- 
citation spectra reflect the presence of a neutral plasma core 
and a positively charged ion shell. 
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Highly excited Rydberg atoms jj) are very sensitive to 
static and AC electric and magnetic fields, and can be 
used for the spectroscopic detection and compensation 
of such fields. Rydberg atoms have been used to detect 
microwave and FIR radiation j^] , to measure and to com- 
pensate weak static electric fields (3) , and to study QED 
level shifts and cavity-QED effects involving only few mi- 
crowave photons . The spectroscopy of Rydberg atoms 
has been employed to diagnose electric fields in DC and 
high-frequency discharge plasmas ||. In this paper, we 
use Rydberg atom spectroscopy to determine the electric 
field in millimeter-sized, strongly coupled plasmas that 
since recently can be generated by the laser excitation 
of laser-cooled clouds of atoms |^] . These plasmas hold 
the promise of bearing novel low-temperature recombi- 
nation Q and oscillation || behavior. Further, there is 
evidence for a Mott transition, by which strongly-coupled 
plasmas and/or dense gases of Rydberg atoms may form 
a new state of highly excited, metastable matter jjj. To 
investigate these and other issues, one could utilize the 
non-invasive field-measurement method employed in the 
present paper. 

Our experiment has recently been described [|l0| . In a 
two-step optical excitation, laser-cooled rubidium atoms 
confined in a volume < 1mm 3 are first excited from the 
ground state 5Si/ 2 to the 5P 3 / 2 state, which has an 
ionization wavelength of Ai on = 479. lnm. A UV laser 
pulse (Atjv = 355nm, duration<10ns, pulse energy up 
to lOmJ) then partially ionizes the 5P 3 / 2 -atoms, produc- 
ing photoelectrons with 0.9eV initial kinetic energy. The 
photoelectrons escape on a time scale of a few ns, leav- 
ing behind a cloud of slow Rb + -ions. About 23ns after 
the UV pulse and while the lower-step laser is still on, 
the remaining 5P 3 / 2 -atoms that are embedded in the ion 
plasma are excited by a blue dye laser pulse (A w 480nm, 
duration<10ns, bandwidth sa 10GHz). The ion cloud 
does not significantly expand before the blue pulse ar- 
rives. The blue pulse excites bound Rydberg states, the 



Stark effect of which we use to measure the plasma elec- 
tric fields. 

We detect the Rydberg atoms using the mechanism 
explained in the following. The ion plasma represents an 
electron trap with a depth Uq given by the ion number, 
the cloud size and the profile of the charge distribution 
||. In most of our experiments Uq < 0.9eV, in which case 
no or only few UV-generated photoelectrons are retained 
in the electron trap (the ones that are moderated to low 
enough energy on their way out). However, the plasma 
electron trap becomes populated with electrons of about 
lOmeV energy due to direct photoionization of 5P 3 / 2 - 
atoms by the ASE with A < Ai on contained in the blue 
laser pulse and due to thermal ionization of a fraction 
of the bound Rydberg atoms excited by the blue pulse 
ifTof . Within a few [is, Z-changing collisions between the 
trapped electrons and the Rydberg atoms efficiently pro- 
mote the latter into long-lived high-(Z) states. During the 
subsequent Coulomb expansion of the plasma, which typ- 
ically takes a few tens of [is, the depth Uq of the plasma 
electron trap approaches zero. During that time, the 
electrons trapped in the plasma are gradually released, 
producing a plasma decay signal in our microchannel- 
plate (MCP) electron detector, which is located about 
10cm from the atomic cloud. A large fraction of the sur- 
viving high-Z Rydberg atoms thermally ionize on a time 
scale of order 10ms |l^] , which exceeds the plasma decay 
time by about a factor of thousand. The delayed ther- 
mal electrons originating in the high-i Rydberg atoms 
are counted in a gate that starts well after the plasma 
decay, yielding a noise-free signal of the Rydberg atom 
excitation. Rydberg excitation spectra are obtained by 
recording the delayed electron counts as a function of the 
wavelength of the blue laser. 

Fig. [l| shows a typical experimental Rydberg excitation 
spectrum of atoms in an electric-field-bearing plasma. 
The field manifests itself through a) the appearance of 
p- and hydrogenic lines, labeled p and h in Fig. [j], b) 
termination of the discrete Rydberg series at a certain 
field-dependent value of the principal quantum number 
n, and c) continuum lowering. While the effects a) and 
b) are solely determined by the state of the ion plasma at 
the time instant of the Rydberg excitation, the feature c) 
also depends on the history of the Rydberg atoms after 
their excitation. Since we intend to analyze the plasma 
at the time instant of the Rydberg atom excitation, we 
disregard the feature c) as a field measurement tool. We 
also don't use b) , because we have found that in many of 
our spectra the laser linewidth terminates the Rydberg 
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series, but not the electric field. 
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FIG. 1. Experimental Rydberg excitation spectrum of 
87 Rb 5P3/2 atoms embedded in a laser-generated ion plasma. 
The three types of spectral ranges with low oscillator strength, 
denoted hs-, mid- and pd-gaps, become filled in at certain crit- 
ical wavelengths that depend on the strength of the plasma 
electric field. The detail plots in the lower row show how we 
determine the critical wavelengths. 

The features of the above type a) develop as follows. 
As n is increased, the electric field first manifests itself 
in the appearance of p-lines and of triangular-shaped 
/i-features (see Fig. [l]), caused by electric- field- induced 
state mixing. Due to the linear Stark effect, the hydro- 
genic states fan out over an energy range 3n 2 E (atomic 
units) ( § and Fig. |). As n increases further, the h- 
features progressively expand and fill in the spectral gaps 
of originally zero oscillator strength between the discrete 
non-hydrogenic s-, p- and d-lines. There are three types 
of such spectral gaps, which we label hs-, mid- and pd- 
gaps. The gaps become filled in with significant oscillator 
strength in exactly that order, at quite well defined crit- 
ical wavelengths. Those wavelengths can be readily con- 
verted into critical effective quantum numbers rths, "-mid, 
and n p d- The rii are robust indicators for the electric 
field, because they are solely determined by the general 
spreading behavior of the quasi-continuous /i-features. 
No absolute line strength or line-strength ratio needs to 
be determined. 

Assuming a constant electric field Eq, the values of 
rij at which the spectral gaps become filled in can be 
related to Eq using the Stark map of Rb. As shown in 
Fig. |^, the /is-gap below the hydrogenic levels with princi- 
pal quantum number n disappears when the energy shift 
of the most red-shifted hydrogenic state, which is approx- 
imately — 3n 2 Eo/2, equals the splitting 0.13n -3 between 
the unshifted hydrogenic energy and the next-lower s- 
state. Further, Fig. [| shows that the mid-gaps, i.e. the 
gaps between the states (n + 3)s and (n + l)d and be- 



tween the field-free hydrogenic states n — 1 and the states 
(n + 2)p, become entirely covered with hydrogenic states 
at practically identical electric fields. Finally, the pd-gaps 
disappear when neighboring hydrogenic manifolds meet, 
i.e. when 2>ti 2 Eq — n~ 3 . The relations between Eq and 
the critical quantum numbers n;n are obtained as 

E = 0.086^ 
E = 0.23n"f do 

Eq = 0.33n p ^ (atomic units). (1) 

The zeros in the subscripts indicate that these relations 
hold for a constant electric field. Atomic species with 
quantum defects different from Rb would follow similar 
relations with different numerical factors. Further, the 
Eqs. share, as expected, the n 5 dependence of the 
familiar Inglis- Teller relation |pd[| . 
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FIG. 2. Stark map of Rb in the vicinity of the n = 19 and 
n — 20 hydrogenic manifolds. 

To model our experimental spectra, we have calculated 
the photoabsorption spectra S(Eq,X) of randomly ori- 
ented atoms in constant electric fields Eq. Since from 
the 5P intermediate state we excite three manifolds of the 
magnetic quantum number |m|, one spectrum S{Eq, A) is 
a properly weighted sum of an m = 0, an m = 1 and an 
m = 2 spectrum. The spectra are obtained as outlined in 
p"2| using a spherical basis set of all states up to n = 90. 
To account for the laser linewidth, the spectra are convo- 
luted with Gaussians of a few GHz FWHM. The calcu- 
lated spectra, one of which is shown in Fig. [| a) , display 
the same qualitative features as our experimental spec- 
tra. This applies, in particular, to the triangular shape 
of the /i-features, which is produced by the unresolved 
excitation of the states in the hydrogenic manifolds. 

An investigation of the calculated photoexcitation 
spectra reveals that the closing of the pd-gap appears 
particularly sudden due to the following peculiarity: for 
parameters between the mid-gap and the pd-gap closing, 
indicated by a double-arrow in Fig. |^, the blue-shifting 
hydrogenic levels within the pd-gap, marked with open 
circles in Fig. carry much less oscillator strength than 
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all the other nearby hydrogenic levels. As a result, the 
whole pd-g&p appears virtually oscillator-strength- free up 
until the pd-g&p closing condition 3n 2 Eo > n~ 3 is satis- 
fied. Once the p<i-gap closing condition is satisfied, the 
oscillator strength of the hydrogenic states inside the pd- 
g&p increases rather abruptly. This behavior enhances 
the clarity of the pd-gap closing. 
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FIG. 3. Calculated isotropic excitation spectra of 87 Rb 5P 
atoms into Rydberg states. Panel a) is for a homogeneous 
electric field Eo = 6 x 10~ 9 at. un.= 31V/cm. Panel b) shows 
the excitation spectrum - upside down - in a Gaussian cloud 
of 4 x 10 5 Rb + ions that has a most probable electric field of 
£?max = 31V/cm. The arrows mark the critical wavelengths 
where the spectral gaps of low oscillator strength become filled 



The electric field experienced by the 5p3/2-atonis is not 
homogeneous but follows a probability distribution P{E). 
Assuming identical Gaussian density profiles for the ions 
and the 5P 3 / 2 -atoms embedded in between them, the dis- 
tribution P(E) can be calculated by a random placement 
of ion charges and 5P3/2-atoms using a radial distribution 
function P(r) = 47rr 2 (27rcr 2 )" 3 / 2 exp(-r 2 /(2cr 2 )). The 
resultant field distributions P{E) have well defined most 
probable electric fields P max , which practically coincide 
with the singularities of the distributions produced by the 
corresponding smoothed charge distributions (compare 
the solid and dotted curves in Fig. ||). The effect of the 
microficlds of the discrete ion charges largely is to smear 
out the singularities produced by the smoothed charge 
distributions, adding fields significantly larger than £ max 
to the field distribution. Deviations of the density pro- 
files of the ions and the 5P3/2-atoms from Gaussians will 
modify P{E), a dominant peak of P{E). However, a 
most probable field P max and a microficld-induced tail of 
fields E > P max will remain. 

We have used simulated electric-field distributions 
P(E) to form averages of the spectra S(Eq, A) weighted 
by P(E). As shown in panel b) of Fig. ||, one can still eas- 
ily obtain critical quantum numbers rihs, "-mid, and n p( j 
at which the respective regions of low oscillator strength 
disappear. The excitation spectra a) and b) in Fig. || 
are, in fact, quite similar, because the P{E) are clearly 
dominated by fields near the most probable one. The 
microfield-induced spread of P{E) from P max towards 



higher fields causes the spectral gaps with low oscillator 
strength to become filled in at values of n that are slightly 
lower than in the case of a homogeneous field E = P max . 
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FIG. 4. Electric-field distributions P(E) in Gaussian ion 
clouds with N = 4 x 10 s and a — 0.2mm. Solid: Distribution 
obtained by random placement of discrete ion charges follow- 
ing a Gaussian density distribution. Dotted: Field distribu- 
tion of the corresponding continuous charge density. Dashed: 
Holtsmark distribution of the microfields in a neutral plasma 
with the same ion density as in the center of the above Gaus- 
sian ion cloud. 

Evaluating a number of spectra as shown in Fig. || 
b), we have determined the relations between the most 
probable fields P max and the critical effective quantum 
numbers nh s , ^mid, and n pc j for Rydberg atoms excited 
in Gaussian ion clouds with a FWHM of about 1mm: 



Pmax = 0.05071!- 5 

P max = 0.145n m f d 

Pmax = 0-26n~j (atomic units), 



(2) 



We have used Eqs. || to experimentally determine the 
values of P max in ion clouds with diameters of ~ 1mm 
and various ion numbers. The spectra obtained in a few 
such experiments are displayed in Fig. ||. The critical 
wavelengths, indicated by arrows, have been determined 
as shown in Figs. [I] and ^, and have been converted 
into critical quantum numbers nhs, n mi d and n p d- Using 
Eq. ^ up to three spectroscopic measurements of P max 
are obtained for each spectrum. The relative uncertainty 
of each single clcctric-ficld value is 30%. Within that un- 
certainty, the different electric-field values obtained for 
the individual spectra agree well. The electric field val- 
ues quoted in Fig. || are the averages over all field values 
obtained for the respective spectra. 

The most probable field P max in Gaussian ion clouds 
can be calculated as P max = (3.1 x 10 _10 A r ion (T^ 2 )Vm, 



with the RMS cloud size a and ion number 7V; C 



We 



have independently deduced -/Vj on from measurements of 
the MOT fluorescence. The values of a could be deter- 
mined from cuts through CCD images of the atom cloud 
fluorescence, yielding values of a ranging from 0.3mm to 
0.5mm. For large atom clouds, we have found that the 
measured fields exceed the calculated ones by up to a 
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factor of ten. As the atom cloud size is decreased, the 
discrepancy drops to about a factor of two. The fields 
deduced from the spectra are, of course, the ones to be 
trusted. The most likely explanation for the discrepancy 
and its dependence on the size of the atomic cloud is that 
the ions were actually not excited in the whole atomic 
cloud, but were concentrated in a significantly smaller 
volume determined by the diameters of the excitation 
laser pulses. 
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FIG. 5. Experimental excitation spectra of 87 Rb 5P3/2 
atoms embedded in photo-excited plasmas that contain the 
indicated numbers iVion of Rb + -ions. The er-parameters of the 
utilized Gaussian clouds of laser-cooled atoms are also shown. 
The most probable electric fields i5 max , shown on the right, 
are spectroscopically determined and carry an uncertainty of 
about 20%. 



shell carries a large electric field. The distribution P(E) 
will thus become bi-modal, featuring a peak near E = 
and a high-field peak. Rydberg atoms excited in the 
shell produce the spectroscopic field signatures discussed 
in this paper, while atoms excited in the core contribute 
a low-field excitation spectrum with resolved Rydberg 
resonances up to large n. The excitation spectra shown 
in the two lowest plots of Fig. |5| appear, in fact, as su- 
perpositions of a high-field and a low-field spectrum, as 
expected for bi-modal distributions P(E). In the low- 
est curve, the high-field component of P(E) apparently 
is strongly dominated by the low-field one, making it 
hard to determine the value of the electric field in the 
shell. The presence of a strong electric field, can, how- 
ever, be guessed from the appearance of p-lines at long 
wavelengths. 

We note that in the case of perfect macroscopic neu- 
trality of the core the microfields in the core would fol- 
low a Holtsmark distribution [[L3| , as displayed in Fig. |J. 
For our charge carrier densities and laser linewidths, the 
Holtsmark fields are too low to be detectable. 

In this work, we have used spectroscopic tools to probe 
the electric fields in small, laser-excited plasmas. The 
technique has a large future potential: using crossed and 
focused Rydberg atom excitation beams, it will be pos- 
sible to perform spatially resolved measurements of the 
plasma electric fields. A variation of the time delay be- 
tween the UV pulse and the blue probe pulse will allow 
one to study the plasma expansion via the decay of the 
electric field. The use of a laser with a narrow band- 
width (~ 1MHz) will allow one to detect smaller fields, 
such as the Holtsmark fields, and to make precise mea- 
surements of the detailed electric- field distribution P{E). 
Measurements of P(E) could reveal inhomogeneities and 
instabilities in the plasma, which are expected to pro- 
duce structures in P(E) at excessive fields. Spatial order 
in the charge carrier arrangement, such as Wigner crys- 
tals H] and the formation of crystalline Rydberg matter 
H may become detectable via measurements of P(E). 

Support by NSF is acknowledged. We thank Prof. P. 
Bucksbaum for inspiring discussions and generous loan- 
ing of equipment. 



If the number of ions in a Gaussian cloud exceeds 
7.8x 10 5 mm _1 cr, the potential depth Uq of the ion plasma 
exceeds 0.9eV, and a fraction of the UV-excited photo- 
electrons will be retained ||. This is clearly the case in 
the lowest two spectra of Fig. [| where we exceed the 
critical ion number by factors > 2.5 and > 10, respec- 
tively. If the fraction of retained electrons is large, a 
macroscopically close-to-neutral plasma core is formed, 
which is surrounded by a positively charged ionic shell. 
The electric field in the core is close to zero, while the 
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